clc;
clear;
A=[3,2,2,2;
   1,2,-2,-2;
  -1,-2,0,-2;
   0,1,3,5];
[v,d]=eig(sym(A))
A1=inv(v)*A*v%这一步是在对原来的矩阵进行对角化
origin=v*d*inv(v)%恢复到原来的函数
v2=[1 0 -1 0;
    -1 0 1 -1;
    -1 -1 1 0;
     1  1 -2 1]
disp('-----------------------')
inv(v2)*A*v2
v2*d*inv(v2)